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ABSTRACT 

Spider is a balloon-borne instrument designed to map the polarization of the cosmic microwave background (CMB) 
with degree-scale resolution over a large fraction of the sky. Spider's main goal is to measure the amplitude of 
primordial gravitational waves through their imprint on the polarization of the CMB if the tensor-to-scalar ratio, 
r, is greater than 0.03. To achieve this goal, instrumental systematic errors must be controlled with unprecedented 
accuracy. Here, we build on previous work to use simulations of Spider observations to examine the impact of several 
systematic effects that have been characterized through testing and modeling of various instrument components. In 
particular, we investigate the impact of the non-ideal spectral response of the half-wave plates, coupling between focal 
plane components and the Earth's magnetic field, and beam mismatches and asymmetries. We also present a model 
of diffuse polarized foreground emission based on a three-dimensional model of the Galactic magnetic field and dust, 
and study the interaction of this foreground emission with our observation strategy and instrumental effects. We find 
that the expected level of foreground and systematic contamination is sufficiently low for Spider to achieve its science 
goals. 

Subject headings: cosmic microwave background, polarization experiments, B-modes, gravity waves, 
analytical methods 



1. INTRODUCTION 

Observations of the cosmic microwave background 
(CMB) have been central to the development of a stan- 
dard cosmological model over the past few decades. A 
key component of this model is inflation, a period of 
accelerated expansion that occurred very early in the 
history of the universe, which enables the model to re- 
produce the flatness and isotropy observed today. Im- 
portantly, inflation also predicts a nearly scale-invariant 
spectrum of primordial density (scalar) perturbations, 
in excellent agreement with recent observ ations of the 
anisotropics in the intensity of the CMB (jLarson et al.l 
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Eonl). 

Although the inflationary hypothesis w as first put for- 
ward several decades ago (|Guthl I1981D , the details of 
the underlying physics that drives the expansion are 
still uncertain. Furthermore, despite the advances made 
through recent observations of the CMB, the constraints 
on the general parameters used to describe infiation are 
weak. As a result, a plethora of plausible inflationary 
scenarios have been advanced, drawing on a wide range 
of proposed physical mechanisms. 

In addition to the density perturbations that seed 
large-scale structure, inflation generates a stochastic 
background of gravitational waves (tensor perturbations) 
that leave a unique imprint on the polarization of the 
CMB by introducing a curl or "i?-mode" component. A 
detection of the _B-mode signature would provide strong 
evidence for inflation and a measure of its energy scale, 
representing a major breakthrough in cosmology. The 
combination of this measurement with our current knowl- 
edge of the scalar perturbations would begin to elucidate 
the dynamics of inflation and place strong constraints on 
the underlying physics. As the current upper limit on the 
energy scale of inflation is around 10^^ GeV, close to the 
GUT scale, measuring the S-mode signature of gravita- 
tional waves in the CMB provides a unique opportunity 
to probe physics at energies far beyond the reach of ter- 
restrial high-energy experiments. 

The amplitude of inflationary gravitational waves is 
parametrized by the tensor-to-scalar ratio, r, which is 
deflned as the ratio of the power in tensor modes to that 
in scalar modes at some pivot co-moving scale, here taken 
to be 0.002 Mpc"^ The combination of CMB tempera- 
ture data with measurements of large-scale structure sets 
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the current upper lim it of r < 0.22 at 95% confidence 
(jKomatsu et al.|[2009D . However, cosmic variance ulti- 
mately limits the r constraints that can be obtained from 
total intensity measurements, and further improvements 
will require direct measurement of CMB polarization. 

A detection of the dominant .E-mode compo- 
nent of CMB po larization was firs t reported around 
eight years ago (jKovac et al.l |2002|) . and since then 
many e xperiments have repor te d further m e asure- 
ments (IReadhea d et al.l 12004 ISievers et al.l 120071: 
Barkats et al. 20051: IBischoff et al.l 120081: ILeitch et al.l 



20051: IMontrov et al.l 120061: llPiacentini et al.l 120061: 
Wu et all 120071: iPage et al.l 120071 ). Recent observations 



have begun to c haracterize detailed features in the E- 
mode spec trum fBrown et alj 120091: iChiang et all 120101: 
IQUIET Co llaboration 2010). However, a measurement 
of the more interesting B-mode component will require 
a substantial improvement beyond present experimental 
sensitivities. 

The weakness of the _B-mode polarization signature 
presents a major experimental challenge. Aside from the 
question of raw sensitivity, there are several other impor- 
tant obstacles that must be overcome. Systematic errors 
must be controlled with unprecedented accuracy if the 
small i3-mode signal is not to be degraded by confusion 
with the much larger E-mode and total intensity sig- 
nals. Diffuse Galactic foregrounds will also contaminate 
the observations, with synchrotron emission and ther- 
mal dust emission expected to significantly complicate 
the interpretation of polarization measurements across 
microwave frequencies . 

Spider (jCrin et al.l [200l iFilippini et all [2010l ) is a 
balloon-borne microwave polarimeter designed to search 
for the imprint of inflationary gravitational waves upon 
the polarization of the CMB. From a vantage point above 
the bulk of the atmosphere. Spider will map the CMB's 
polarization over a large portion of the sky with a res- 
olution close to one degree. Spider's main goal is to 
measure the amplitude of gravitational waves with sen- 
sitivity down to r = 0.03 (with 99% confidence), or to 
place an upper limit at this value if no detection is made. 
To obtain this sensitivity. Spider will employ over 2500 
b olometers. 

iMacTavish et al.l (|2008D investigated the potential im- 
pact of many common sources of systematic error on the 
ability of Spider to achieve its main scientific goal, set- 
ting tolerance limits to guide the design of the instrument 
and observation strategy. In this paper, we continue this 
work, moving our focus onto new systematic effects that 
have proven to be important but that have so far re- 
ceived less attention. Instead of setting tolerance limits, 
our aim is to ensure that the measured performance of 
various instrument components and prototypes is suffi- 
cient to meet our main goal. 

We also develop a detailed model of the foreground 
emission Spider will see. This makes the sky model 
against which the systematic effects are assessed more 
complete, and also enables us to further test the robust- 
ness of our observation and analysis strategy. In this 
work we do not attempt to address the important issue 
of separating the foreground and CMB contributions to 
our data. 

This paper is organized as follows. Section [2] gives 
a brief overview of the instrument and introduces the 



systematic effects we investigate. The model of polar- 
ized foregrounds we use is described in Section [3l Sec- 
tions |4] and [5] introduce the simulation methodology and 
the baseline observation strategy we use, and present the 
results of these studies. Finally, Section [6] gives a discus- 
sion of the results and our conclusions. 

2. THE INSTRUMENT 

A detai led description of the Spider instrument is pre- 
sented in iRunvan et al.l ()2010[ ). Here we provide a brief 
overview of the instrument and the sources of error we 
consider. 

The Spider payload consists of six separate monochro- 
matic instruments housed in a single liquid-helium cryo- 
stat. All six point in the same direction on the sky. The 
optical de sign is based on that o f the successful Bicep 
telescope (jTakahashi et al.l 12010"). Each instrument in- 
cludes a telecentric refracting telescope consisting of two 
lenses cooled to 4K. At 150 GHz this design produces 
a beam on the sky with full-width at half-maximum 
(fwhm) close to 40arcmin. Each telescope focuses light 
onto a focal plane array of antenna-coupled bolometers, 
which are cooled to 300 mK using '^He closed-cycle sorp- 
tion fridges. Each dual-polarization spatial pixel in the 
focal plane consists of two phased arrays of slot dipole 
antennas, which are sensitive to polarization along or- 
thogonal axes. The power incident on these antennas is 
deposited on to a bolometer and detected with a super- 
conducting transition-edge sensor (TES). Superconduct- 
ing quantum interference devices (SQUIDS) are used to 
amplify the current produced by the TES s, which is then 
read out using time-domain m ultiplexing ()de Korte et al.l 
120031: iBattistelh et al.l [20081 ). The fo cal-plane a r chitec - 
ture is described in more detail in iKuo et al.l (|2008[ ). 
Note that the design of the beam-defining optical com- 
ponents, the lenses and the focal plane a.rrays, is very 
simil ar to that of the Bicep2 experiment (jOgburn et al.l 
l2010f ). although the ar chite cture of the focal plane unit 
is different; see Section [5731 

A rotatable cold half-wave plate is mounted at the 
aperture of each telescope, on the sky-side of the primary 
lens. The plates are constructed from 330-mm diam- 
eter slabs of birefringent single-crystal sapphire, with a 
quarter-wave quartz coating to minimize reflections. The 
orientation of the half- wave plate determines the orienta- 
tion of the instrument's polarization sensitivity axes with 
respect to the sky (see equation ([6])). By making observa- 
tions with the half- wave plate oriented over a range of an- 
gles during the flight, high polarization cross- linking can 
be achieved, improving the fldelity of the recovered polar- 
ization maps. The half-wave plates also provide powerful 
defense against several important sources of systematic 
error. Any error that does not transform as a true polar- 
ization signal under rotation of the half-wave plate can 
be distinguished from CMB signal during data analy- 
sis. This suppression will occur naturally using standard 
map-making techniques (provided that the polarization 
cross- linking is sufficient), as the systematic contribu- 
tions to the data will have functional forms nearly or- 
thogonal to that of the true sky signal. If necessary, 
further suppression may be achievable by modifying the 
map-making algorithm to model the contribution to the 
data induced by the error, and marginalizing over it if 
required. 



Spider Optimization II 



3 



In this paper, we consider several potential sources of 
systematic error in the Spider instrument. First, the be- 
havior of the half-wave plate is weakly dependent on the 
frequency of the incident radiation, which leads to non- 
idealities when averaged over the Spider bandwidths, 
which are around 25% of the central frequency. Second, 
the focal plane is highly sensitive to external magnetic 
fields, and the motion of the instrument through the 
magnetic field of the Earth has the potential to generate 
significant spurious signals. Using laboratory measure- 
ments of the magnetic pickup in the focal plane, we esti- 
mate the in-flight response to the Earth's field. Finally, 
we consider imperfections in the optical beams defined 
by the antenna arrays and lenses, focusing on asymme- 
tries and mismatches between the orthogonally polarized 
detectors in each focal plane pixel. Optical characteriza- 
tion of the Spider and Bicep2 instruments have shown 
that these beam systematics are present at a small level, 
and we use these measurements to guide our studies of 
the expected impact on Spider's science goals. 

3. SKY MODEL: FOREGROUNDS 

Galactic foreground emission is expected to contribute 
significantly to the polarized microwave emission across 
the sky, and may well dominate the CMB gravitational- 
wave signal at all frequencies. Therefore, it is impor- 
tant to include the foreground emission in our sky model 
for this work. Furthermore, Spider will observe a large 
fraction of the sky, potentially including regions close to 
the Galactic plane where the foreground emission is very 
bright in comparison to the CMB. The presence of this 
bright emission in our data may affect the performance 
of the observation and analysis strategy we use. 

Unfortunately, these foregrounds are poorly con- 
strained by current data and poorly understood. This is 
particularly true above around 90 GHz, where the CMB 
emission is strongest and where Spider will operate. At 
these frequencies, the foreground emission is expected 
to be dominated by thermal emission from interstellar 
dust. Direct measurements by Archeops (Benoit ct all 
2004:'Pont hieu et al.ll2005f ) and by WMAP (Kogut_eLd] 
2007) have shown that this emission is polarized at the 
few percent level at frequencies higher than 90 GHz. 
This is in agreement with our theoreti cal understanding 
of th e underlying physical mechanism (jPraine &: Fraissg 
[2009h . which' can be summarized as follows. Dust grains 
are generally non-spherical, and preferentially emit ra- 
diation polarized along their longest axis. Mechanisms 
exist that align these grains with this axis perpendicular 
to the Galactic magnetic field. This leads to a net lin- 
ear polarization of the observed emission. Independent 
observations of the polarization of starlight (Heilcs 199^; 
iFosalba et al.l[200 2!) are consistent with this picture. 

Synchrotron emission, generated by the gyration of 
cosmic ray electrons in the Galactic magnetic field, is in- 
trinsically polarized and constitutes the main polarized 
foreground at lower frequencies (Page et al. 2007). Since 
the polarized emission from interstellar dust is expected 
to dominate the polarized Galactic synchrotron emission 
at frequencies above ^ 70 GHz, we do not include the 
latter in our sky model. 

Although little data is available regarding the polarized 
emission from dust, the same is not true of its total inten- 
sity. In particular, the IRAS satellite observed this emis- 



sion across the sky at 100 /im and 240 /im, close to the 
peak in the dust emission. By constraining physically- 
motivated extrapol ations of t hese o bservations using fur- 
ther data, Finkbei ner et al.l (|1999l hereafter FDS) pro- 
vided models of the emission at microwave wavelengths. 
At 94 GHz, these models have been shown to agree well 
with the WMAP observations, up to a small normal- 
ization factor and s ome minor struct ural differences in 
the Galactic plane ()Gold et al.ll2009[) . Data are more 
limited in the higher frequency bands relevant to Spi- 
der, bu t measurements agree well with the FDS pre - 
dictions (jCulverhouse et al.ll2010trV'eneziani et al.ll2010H . 
We use FDS's model 8 to trace the total intensity of the 
dust emission. 

The other necessary fields for a full description of the 
dust emission, namely its degree and direction of polar- 
ization, are highly dependent on the Galactic magnetic 
field. Since the observed polarization signal is the sum 
of the emission from many independent regions along the 
line-of-sight, it is sensitive to the three-dimensional struc- 
ture of this magnetic field. As a result, modeling the de- 
gree and direction of polarization of the dust emission re- 
quires the evaluation of appropriate linc-of-sight integrals 
given three-dimensional models of the Galactic magnetic 
field and of the other relevant Galactic constituents. 

Away from the Galactic center, the Galactic magnetic 
field is usually considered to have two nearly independent 
components: a large-scale coherent field associated with 
the Galactic disk and a small-scale field arising from tur- 
bulence in the interstellar plasma sourced by astrophys- 
ical events such as supernovae and stellar winds. The 
most informative probes of these fields are Faraday rota- 
tion measures of p ulsars and extra-Galactic radio sources 
(iHavcrkor n et al.l 120061: iHan et al.l I2006D . Whilst there 
is general agreement that the large-scale field follows a 
spiral pattern, its detailed structure is still uncertain. 
This uncertainty is unimportant when considering areas 
of sky at high Galactic latitudes: since the dust is con- 
centrated in a thin disk about the Galactic plane, we 
only see emission originating within around 1 kpc or so 
of the Sun, a region in which the large-scale field is rea- 
sonably well characterized. However, as Spider's large 
sky coverage may include part of the Galactic plane, we 
require a model of the large-scale field structure in that 
plane. One popular candidate is th e bisymmetric spiral 
i|Han fc Qiaolll994l : iSun et al.l [20081 ) . which can be writ- 
ten as 

Bp = —Bq cos { (j) + ipln—] sinp cos yq 

V PoJ 

Bj, — —Bq cos ( -f -0 In — I cosp cos xo 

V PoJ 

B^ = Bosmxo- (1) 

Here, p, (jj and z are Galactocentric cylindrical coordi- 
nates with (j) measured from the direction of the Sun, 
p is the pitch angle of the field, ^ = 1/tanp, po de- 
fines the radial scale of the spiral, and xo parametrizes 
the amplitude of the z component. We use the param- 
eter values suggested in iMiville-Deschenes et al.l ()2008D : 
p = —8.5 degrees, po = Hkpc and xo = 8 degrees, with 
the field amplitude set to So = 3 pG and the distance be- 
tween the Sun and the Galactic center taken to be 8 kpc. 
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Although this field model i s unlikely to provide a full 
descr iption of our Galaxy Men et alJ [20081: [S un et alj 
it is sufficient for our current purpose since we 
do not require an accurate template of the sky, only a 
reasonable approximation against which to test the per- 
formance of the experiment. 

The turbulent component of the Galactic magnetic 
field is somewhat less well unde rstood. When constrain- 
ing th e above large-scale field, iMiville-Deschenes et al.l 
(|2008| ) simultaneously fit a small-scale field component 
with best-fit r.m.s amplitude -Br.m.s. = 1-7 fJ-G. Since 
several different studies agree that the r.m.s. amplitude 
is similar to the amplitude of the large-scale field i n 
the solar vicinity (iFosalba et all 120021: iHan et al.ll2006D . 
we set Sr.m.s. = 2 nG. iMinter fc Spangled (llQQel Tex- 
amined the rotation measures of extra-Galactic sources 
across a small patch of sky and concluded that the 
data were consistent with Kolmogorov turbulence on 
scales smaller than 4pc, assuming a statistically isotropic 
and homogeneous Gaussian field. On larger scales, 
they found a somewhat flatter energy spectrum with 
an outer scale of up to 96 pc. Kolmogorov-type spec- 
tra up to kilo-parsec scales in the interstellar mag- 
netic field and other interstellar plasma components have 
also been reported by other studies (lArmstrong et al.l 
[TOOllLazarian fc Pogosvanl2000HCho fc Lazarianll2008[ ). 
Since it is numerically intractable to generate a full-sky 
realization of this turbulent field in three dimensions at 
sufficiently high resolution, we resort to independent one- 
dimensional realizations along the line-of-sight to each 
pixel. This model ignores correlations across the sky, 
but properly incorporates the line-of-sight depolariza- 
tion. We choose an injection scale of 100 pc, assume the 
dissipation scale to be small, and use the one-dimensional 
Kolmogorov energy spectral index of —5/3. 

We model the large-scale spatial distribution of the 
dust densit y, Ud, using a simplification of the model con- 
strained in lDrimmel fc Spergell ()200lD . 



nd = riQ exp 



Pd 



sech" 



(2) 



For consistency w ith the WMAP polarization analysis 
(jPage et al.ll2007[) . we take the scale height Zd — 100 pc 
and the scale radius pd — i kpc. We do not attempt 
to model the small-scale variations in the dust density 
and temperature, which may also affect the polarization 
degree and direction. Small-scale variations in the total 
intensity are included via the FDS model. 

The model also requires a description of the physics 
of grain alignment and of the intrinsic polarization of 
the emission from an individual grain. In general, these 
are complex functions of the magnetic field and of var- 
ious properties of the grains. Recently, good progress 
has been made in describing the deta ils of the alignment 
using the theory of radiat ive torques (jLazarian fc Hoan"^ 
l2007HHoang fc Lazaria3l2008[ ). However, it is still diffi- 
cult to produce a well con strained qu antitative descrip- 
tion to apply to our model (jLazarian fc H oang 2009). We 
instead describe the alignment in an integrated manner, 
without recourse to the details of a particular physical 
mechanism. We assume that the polarization direction 
is always perpendicular to the component of the mag- 
netic field in the plane of the sky, and that the degree 



of polarization depends quadratically on the magnetic 
fie ld strength . This is similar to the behavior assumed 
in iPage et al.l (f2007t ). 

We compute the Stokes parameter maps associated 
with our three-dimensional model using the appropriate 
line-of-sight integrals, 

/3d(0,0) =eM / "'^^dMdr 



Jo 

Q3b{0. 0) = I nd{r)pQ[B^{rf - Be{rf] dr 
Jo 

Um{e. 4>) = [ nd{r)po[2B4r)Be{r)] dr, (3) 
Jo 

where the normalization pq is chosen to reproduce 
the average polarization fraction o f 3.6% reported b y 
WMAP outside their P06 mask (jKogut et all I2007D . 
Here, r, 9 and (j) represent spherical polar coordinates 
centered on the telescope, and e is the emissivity of the 
dust as a function of frequency, v. Note that we con- 
form to t he default convent ion applied in the healpix[3 
package (jGorski et al.ll2Q05l) regarding the sign of U. 

From this model we require maps of the polarization 
direction, 7, and degree, P, which are given by 



P{0,, 



(Q 



2 

3D 



71 



J 3D 

1 , fUw 
— - arctan — — 

2 \Qw 



(4) 



Finally, we scale by the FDS intensity maps to obtain 
our final dust model expressions: 



Qdustl^*, 0) = /fds^cos27 
Udust{0,(l)) = /FDS^'sin27. 



(5) 



Figure [T] shows the Stokes parameter fields for this 
model evaluated at 150 GHz. In order to compare the 
dust emission to the CMB, we plot the pseudo-power 
spectra of these maps (calculated using the healpix 
anafast facility) in Figure [21 masking regions within 
10 degrees of the Galactic plane. 

4. SIMULATION METHOD 

We use the simulation pipeline and nominal observa- 
tio n and analysis scheme introduced and investigated 
in IMacTavish et al.l (j2008D . This in turn is based 
largely on the analysis pipeline develope d for the 2003 
flight of Boomer ang |Masi et aTl 120061 ). described in 
iJones et al.l ()2007[ ). 

To limit the computational requirements, data are sim- 
ulated for only 8 spatial focal-plane pixels (16 bolome- 
ters) in a single, evenly-spaced column spanning one Spi- 
der focal plane. The two antennas in each spatial pixel 
are sensitive to polarization parallel and perpendicular 
to the column. This limited simulation is sufficient for 
our current purpose, where we perform signal-only simu- 
lations and investigate systematics effects which are un- 
likely to worsen as the number of bolometers considered 
increases. Therefore, the contamination generated by the 

^* See |http://healpix.jpl.nasa.gov| 
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Fig. 1. — Stokes parameter maps (from top /, Q and (7) in Galac- 
tic coordinates for our model of thermal dust emission at 150 GHz. 
The grey lines show the region observed with the observation strat- 
egy described in Section |4] which excludes the Galactic center. 



various systematic errors reported here can be treated as 
reasonable upper limits on the likely residuals in the ex- 
periment itself. 

We simulate a mid-November launch from Alice 
Springs, Australia (longitude 128.5 degrees east, latitude 
25.5 degrees south), with the balloon subsequently drift- 
ing in longitude at 3.76 x 10~^ degrees per second. The 
gondola is spun at 36 degrees per second at a fixed ele- 
vation of 49 degrees, and we observe only when the Sun 
is ten degrees or more below the horizon. The half-wave 
plate is stepped once per day by 22.5 degrees, and so we 
consider a four day observation period which allows us to 
fully sample the polarization modulation the wave plate 
induces; see equation ([6]). With this scanning strategy, 
Spider observes approximately 60% of the sky. 

Based upon this scenario, the flight simulator gener- 
ates the time-ordered right ascension and declination ob- 
served by each spatial pixel, as well as its orientation 
angle '0 with respect to the sky basis. The sky model is 



generated using the healpix synfast facility and the fore- 
ground model described in Section [3] using nsidc = 1024, 
i.e. pixels of size 3.4arcmin. The sky is also smoothed 
with a Gaussian kernel with fwhm of 40 arcmin to ac- 
count for the angular resolution of Spider. The model 
is generated at 150 GHz, Spider's main science chan- 
nel. As we are not concerned with the detailed impact of 
instrumental noise or foreground separation, this is suf- 
ficient for this work. Note that we use the same CMB 
realization for each of the various simulations we per- 
form. The pointing information is then used to generate 
the time-ordered data from the sky model and the wave- 
plate angle, 0. The bolometer outputs are given by 

d± —Isky ± [Qsky COS 2(V' + 20) + C/sky sin2(?/; + 20)], 

(6) 

where the signs refer to the different bolometers in a 
spatial pixel. Including systematic effects requires mod- 
ifications to the above procedures, which we describe in 
Section [5] for the systematic errors we consider. 

Spider maps are then made using an adap tation of the 
Jacobi iterative map-maker as described in iJones et al.l 
(j2007t ). Although the simulations are noiseless, the maps 
are made using a noise covariance matrix based on a 1// 
spectrum with knee frequ ency lOOmHz a typ ical value 
for the Spider detectors (|Runvan et all 120101 ). We also 
apply a one-pole high-pass filter at 10 mHz to the time- 
ordered data during map-making, as this is expected to 
be necessary for the real data to mitigate long-time scale 
systematic errors such as incomplete knowledge of the 
system transfer functions. Thus we properly include the 
expected loss of modes due to the low frequency noise and 
filtering and the resulting degradation of the recovered 
maps. These maps are produced at ngido — 256, i.e. a 
pixel size 13.7 arcmin. 

By examining maps of the residuals, i.e. the difference 
between the output Spider maps and the input model 
sky, we can begin to evaluate the consequences of the 
various sources of error. In order to properly assess their 
impact on our ability to meet our main goal of measuring 
the gravitational- wave signal, however, it is necessary to 
proceed to the power spectrum domain, where we can 
make direct comparisons with the expected CMB spec- 
tra. 

We estimate spectra from our maps via a spherical har- 
monic transform using anafast. The decomposition of po- 
larization fields into E- and _B-modes is not unique on the 
cut sky, however, and so these so-called 'pseudo'-spectra 
estimates contain mixing between the two. The result 
of this is that the recovered i?-mode spectrum is biased 
by the much larger inpu t i5-mode spectrum (Bunn et aTl 
!I003; Lewi s et al.l [2001 . 

If unaccounted for, this E-B mixing dominates the 
systematic contribution to the recovered i?-mode spec- 
trum. The development and implementation of an op- 
timal unbiased power spectrum estimator which prop- 
erly accounts for this mixing is beyond the scope of this 
work. Instead, we eliminate this bias by estimating the 
pseudo-spectra of the residual maps. This removes the 
true sky from our estimates, allowing us to concentrate 
on the contribution from systematic errors. We also lose 
a small contribution to the recovered _B-mode spectrum 
from map-level errors which correlate with the sky (e.g. 
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Fig. 2. — Power spectra (from left I, E and B) of our model of thermal dust emission at 150 GHz. The CMB spectra for a universe 
consistent with current data with r = 0.03 are shown for comparison, (red, long-dashed). In the right-hand panel, the contribution from 
gravitational Icnsing is also shown (green, dashed). 



an absolute calibration error), but fully retain the much 
more worrisome systematic effects which couple the re- 
covered i?-modes to i?-modes and total intensity. Using 
the power spectrum of the residuals has the added ad- 
vantage of eliminating the cosmic variance contribution 
from the sky signal, allowing the systematic-induced bias 
to be accurately determined without the need for many 
independent simulations. 

The resulting power spectrum estimate, Cj^''^"^, should 
be compared to the expected CMB B-mode pseudo- 
spectrum measured by Spider without foregrounds, sys- 
tematic errors or E-B mixing. We estimate this by 
performing a systematics- and dust-free simulation with 
no input i?-modes and calculating the resulting pseudo- 
spectrum, C;^'"°^. Finally, to aid interpretation we form 



Mask: Normalized observations per pixel 
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C 



B,noE 



(7) 



so that the final comparison can be made against the 
input CMB spectrum, Cf . 

We emphasize that this procedure is not designed to 
give a complete picture of Spider's ability to estimate 
the i3-mode power spectrum but to allow us to ensure 
that the systematic-induced bias is small compared to 
the expected signal of interest. 

For the above power spectrum estimations we use a 
mask covering the region observed, weighted by the time 
spent observing each sky-pixel to reduce the impact of 
poorly sampled pixels near the map edges. We also add 
a ±10 degree Galactic-plane mask, as this area is heavily 
obscured by foregrounds and will not be available for 
cosmological analysis. This mask is shown in Figure |3l 

5. SIMULATION RESULTS 

5.1. Dust model 

We begin by investigating whether the addition of dust 
to our sky model reduces the effectiveness of the nom- 
inal observation and analysis strategy described in Sec- 
tion m To do so we perform two simulations, one using 




Fig. 3. — Map showing the sky-mask used when calculating power 
spectra in our simulation pipeline in equatorial coordinates. 



our full sky model and one with a CMB-only sky. Figure 
|4] compares the residual sp ectra, Ri, from these sim ula- 
tions. As demonstrated in iMacTavish et aLl ()2008D . ten 
map-making iterations is sufficient to reduce the resid- 
uals to a negligible level for the CMB-only case. Us- 
ing our updated sky model, the residuals remain signif- 
icant after ten map-making iterations, with Ri greater 
than the CMB S-mode spectrum for r — 0.03 on large 
scales {I < 20). We find that increasing the number 
of iterations to around 30 suppresses the residuals to a 
negligible level. Figure [S] shows the residual Q map for 
this case (the residual U map has similar amplitude), 
and demonstrates that the recovered maps are of high 
quality, with little loss of information arising from the 
time-stream filtering and 1// noise. Note that the small 
Galactic regions with relatively high residuals are behind 
our Galactic mask. In accordance with these results, we 
use 30 map-making iterations for the remaining simula- 
tions reported in this paper, unless otherwise stated. For 
completeness, Figure S] also shows Ri after 60 iterations, 
at which point the map-maker has fully converged. 

5.2. Non-ideal half-wave plate 

With the sky model thus established, we turn to our 
investigation of the impact of instrumental systematics 
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Fig. 4. — Residual _B-mode spectra, i?;, for simulations run with and without the dust emission (green and red respectively) after ten 
iterations of the map-maker. Residuals for simulations with dust emission after 30 (dashed) and 60 (long-dashed) iterations are also shown. 
The CMB _B-mode spectrum for r = 0.03 is shown for comparison (black). This confirms that our previous nominal observation strategy 
is not compromised by the inclusion of the dust emission, with residuals orders-of-magnitude below the level of the signal of interest after 
30 iterations. 



Residual Q: nominal case (with dust) 




Fig. 5. — Residual Q polarization map for a simulation run with our complete sky model (CMB and dust) after 30 iterations of the 
map-maker in equatorial coordinates. The black lines show the Galactic cut used when calculating power spectra. The residuals are low 
across the observed sky apart from some small regions in the Galactic plane. 



on Spider's expected sensitivity. We begin with the non- 
trivial spectral response of the Spider half-wave plate. 

The Spider half-wave plates consist of birefringent 
sapphire with thickness chosen such that the optical path 
length difference along the crystal axes is exactly a half 
wave-length at the center of the appropriate frequency 
band. The anti-reflection quartz coating is also opti- 
mized to this frequency. However, as we move away from 
this central frequency, the plate behavior becomes non- 
ideal; the optical path difference is no longer exactly a 
half wave-length, and the amplitude of reflections, al- 
though still suppressed, increases. When averaged over 
the frequency band these non-idealities remain. Further- 
more, the band-averaged non-idealities depend on the 



emission spectrum, S{iy), of the incident radiation. This 
spectrum is a priori unknown, as we do not know the rel- 
ative contributions from the CMB and the various fore- 
grounds. This complicates efforts to predict the effects 
of the half-wave plate non-idealities, and to account for 
them during data analysis should it prove n ecessary. 

We u se the analytic formula derived in iBrvan et al.l 
(|2010bf ) to calculate the band-averaged Miiller matrix de- 
scribing the Spider half-wave plate for both CMB and 
dust emission spectra. This formula is based on a phys- 
ical optics model of the half-wave plate. We use the 
refractive index measurements reported in IBrvan et al.l 
(2010a) for sapphire at 5K, and take the quartz anti- 
reflection coatings to have a refractive index of 1.95. In 
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T 


P 


c 


s 


Ideal 


1 





-1 





CMB 


0.97127 


0.00673 


-0.94157 


0.03126 


Dust 


0.96952 


0.00628 


-0.94067 


-0.01190 



TABLE 1 

Parameters describing the band-averaged half-wave plate for a 
Spider detector response function, and the expected emission 
spectra for the CMB anisotropics and thermal dust. Values are 
also given for an ideal half-wave plate. 



performing the band-averaging, the dust emission spec- 
trum is modeled foUowing model number eight in FDS, 
using the average temperatures for the two model dust 
components, Ti = 9.7K and T2 — 16.2K, and so 



'S'cMB(i^) — 



•S'dustCl') 



dB{T, j^) 



dT 



TcME 



1^0 



1.67 



S(Ti,iy) + 0.933 



2.7 



S(T2,J/), 

(8) 



where B{T, v) is the black-body spectrum and is 
150 GHz. We also incorpor ate the Spider detect or re- 
sponse spectrum reported in lFilippini et al.l (|2010l ). The 
band-averaged half-wave plate Miiller matrices are fully 
described by four independent parameters such that 



M 



HWP,band 



(T 

P 




(9) 



Some of these parameters have straight forward interpre- 
tations. T represents an overall transmission loss, and p 
arises due to differences between the transmission spectra 
for radiation polarized along the two crystal axes. The 
calculated parameter values for our 150 GHz band are 
given in Table [TJ although note that as the CMB is ex- 
pected to have no circular polarization, the s parameter 
should not be relevant to Spider. 

The Miiller matrix describing the transmission of ra- 
diation through the entire instrument from the sky to a 
bolometer can be expressed in terms of Mnwp.band and 
the matrix describing the antenna polarization, Mant, as 



Mi 



=M«ntR.R-flM 



6/ IVlHWP,band Re Ri/. , 



(10) 



where e is the angle between the antenna's polarization 
direction and ip i.e. its orientation in the focal plane (here 
taken to be and 90 degrees for the two bolometers in a 
spatial pixel), and the matrix Rq, describes a basis rota- 
tion through an angle a. We assume the polarization of 
the antenna is perfect and so 



Ma 



f 1 1 \ 

110 
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(11) 



The bolometers measure the total power incident upon 
them, and so the ideal bolometer outputs, equation (|6]), 



are replaced by (|Brvan et al.ll2010bl) 

±(T - c)[Qsk.v cos2(i/; + 20) + C/sky sin2(7/; + 261)] 
±2p[/sky cos 29 - Qsky cos 2(V' + 0) - C/^ky sm2{ilj + 6)] 
±{T + c)[Qsky cos2?A + Usky sin2V']. (12) 

To add the half-wave plate non-idealities to our simu- 
lations we make this replacement when generating the 
time-ordered data. 

Examining equation (jl2[) . we see that the total in- 
tensity contribution is reduced by the transmission loss 
caused by the half-wave plate, as expected. Similarly, 
the polarization efficiency is reduced to (T — c)/2. Both 
of these effects are degenerate with the more general cal- 
ibration requirements of the experiment. Of greater in- 
terest are the remaining contributions to the bolometer 
outputs, which have dependencies on the half-wave plate 
orientation unlike those of the ideal outputs, being either 
independent of 6 or sinusoidal in 29. As these new contri- 
butions transform differently under rotation of the half- 
wave plate to the ideal bolometer outputs, these errors 
should be suppressed at the map level by the half-wave 
plate modulation. However, our nominal four day scan 
is not sufficient to fully sample these terms to take full 
advantage of this suppression, and so for this section we 
increase our scan duration to eight days. 

In the following simulation results, we calibrate the 
Spider maps produced by the pipeline by applying the 
appropriate factor [for polarization 2/(T — c)] before cal- 
culating the residual maps and spectra. This enables 
us to investigate the more complicated and less well un- 
derstood systematics introduced by the half-wave plate. 
Although the calibration factor varies between the sky 
components, the difference is small (0.1%) in compari- 
son to our absolute calibration requirements and to the 
likely level of foreground residuals after foreground sep- 
aration. 

As the half-wave plate parameters vary between the 
sky components, we perform two independent simula- 
tions for the CMB and dust. The residual spectra for 
these simulations are shown in Figure IHl For the dust- 
only simulation the residuals are not negligible, rising to 
over 40% of the expected CMB B-mode spectrum for 
r = 0.03 at low /. The CMB-only residual is lower, but 
rises to around 20% of the expected CMB spectrum. The 
residual Q map for the dust-only simulation is shown in 
Figure [7] (the residual U map has similar amplitude). By 
running a further dust-only simulation on an unpolar- 
ized sky (i.e. with Q = U = everywhere) and setting 
T and c to their ideal values, we see that the residuals 
are dominated by leakage of the dust total intensity into 
polarization via p, as shown in Figure [51 Similar simula- 
tions isolating the contribution from the {T + c)/2 terms 
in equation (|12p confirm that such an error does not lead 
to significant residuals for any reasonable parameter val- 
ues. Note that after 30 iterations, the map-maker has 
fully converged for these simulations. 

Note that taking the unrealistic step of masking the 
Galactic plane on the input sky significantly reduces 
the residuals resulting from the half-wave plate non- 
idealities. Although the plane is always masked in the 
Spider maps when we calculate spectra, the filtering of 
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Fig. 6. — Residual _B-mode spectra, Ri, for simulations using the non-ideal half-wave plate model for the CMB (magenta, dashed) and 
dust (green, long-dashed). The spectrum for our nominal simulation (red) and the CMB B-mode spectrum for r = 0.03 (black) are shown 
for comparison. The spectrum for a simulation with p = Pdust E'ud only the dust total intensity on the sky (i.e. no polarization) is also 
shown (dotted blue). 



Residual Q: HWP model for dust emission 




Fig. 7. — Residual Q map for a dust-only simulation including the half-wave plate model in equatorial coordinates. The black lines show 
the Galactic cut used when calculating power spectra. 



the time-ordered data partially delocalizes the large sys- 
tematic signal generated as the instrument scans across 
the plane, moving the resulting residuals to high latitude 
regions of the map. 

Our simulations have shown that the errors introduced 
by our non-ideal half-wave plate model have the poten- 
tial to mildly compromise Spider's science results if not 
accounted for during map-making. However, the relevant 
parameters describing the half-wave plate non-idealities 
are not a strong function of the emission spectrum of the 
sky: comparing the values for the CMB and average dust 
spectra, the most troublesome parameter, p, changes by 
only around ten percent, and T and c change by less 
than one percent. Similarly, variation in these parame- 
ters across the sky (as, for example, the dust tempera- 
tures vary) are expected to be smaller still. Therefore, 



we expect to be able to accurately correct for these er- 
rors during map-making by upgrading the algorithm to 
include the extra contributions to the bolometer outputs 
in the pointing matrix, using constant cal i brated values 
of T, p and c, as proposed in iBrvan et al.l (|2010bD . This 
should reduce the map-level residuals by at least an order 
of magnitude, i.e. to a negligible level given our current 
science goals. For future experiments targeting r < 0.01, 
more detailed simulations will need to be undertaken to 
investigate the viability of this strategy. Modifying the 
map-maker in this manner may also provide greater flexi- 
bility in Spider's scanning strategy, as the increase from 
a 4-day cycle in our wave plate angle to an 8-day cycle 
might not be required. 

5.3. Magnetic field sensitivity 
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A second potential concern for Spider is stray mag- 
netic field pickup in the superconducting electronics 
chain. The SQUID amplifiers used in the signal read- 
out chain are essentially highly sensitive magnetometers; 
left unshielded, they will not only amplify signal from 
the bolometers but also from the Earth's magnetic field. 
Furthermore, the transition temperature, Tc, of the TES 
detectors is weakly dependent on the magnetic field en- 
vironment. Motion of the gondola through the terres- 
trial magnetic field may thus induce spurious signals in 
the Spider time streams. In order to mitigate these ef- 
fects, considerable effort has been made to shield the 
focal plane so that pickup fr om the Earth' s field will 
be sub-dominant to the CMB. iRunvan et al.l (|2010[ ) de- 
scribes recent changes to the Spider design, notably a 
move from a flat focal plane architecture to a shielded 
box scheme, with the aim of significantly reducing the 
magnetic pickup. 

These improvements have been successful, and the 
magnetic response of a Spider telescope is now low 
enough that it has proven difficult to characterize 
through laboratory measurements. Response amplitudes 
have been measured for many of the detectors for the 
three magnetic field axes, with the remaining amplitudes 
being too small to detect against the measurement noise. 
These measurements were taken with the TES bolome- 
ters on transition (as they would be in flight), so that any 
changes in either Tc or SQUID fiux would be detected. In 
order to assign amplitudes to the eight detector pairs in 
our simulation, we have selected randomly from among 
the detected response amplitudes, taking care to ensure 
that the variation in amplitudes across the focal plane 
and across the different magnetic field axes is typical of 
that seen in the measurements. 

It is not only the overall amplitude of the pickup that 
is important; any difference in magnetic response within 
a detector pair will lead to spurious polarization signals. 
We have therefore preserved detector pairings when se- 
lecting the amplitudes for our model. However, as these 
differences are typically much smaller than the overall 
amplitude, they are not well constrained by our mea- 
surements. In reality the differences may be significantly 
smaller, and so the simulations presented here represent 
an upper limit on the spurious signal Spider is likely to 
see. 

Due to the very low signal levels, characterizing the 
frequency dependence of the pickup is not feasible, so 
instead we use a spectrum typical of those measured for 
the previous focal plane design: a one-pole filter with a 
time constant of 0.3 s. 

To include magnetic pickup in our simulations, the re- 
sponse amplitudes and modeled frequency dependence 
are used to convert the input Earth's magnetic field (in 
Tesla) to the output stray pickup observed in the Spider 
time streams (in /iK). The Earth's magnetic field is deter- 
mined for each Spider pointing (latitude, longitude, alti- 
tude and time) using the World Magnetic Model code and 
data provided by the National Geophysical Data Cen- 
teP^. which estimates the strength and direction of the 
Earth's main magnetic field for a given point. We use the 
pointing time streams to account for the changing orien- 
tation of the focal plane with respect to the terrestrial 

http:/ /www. ngdc.noaa.gov/geomag/WMM/ 



field. In order to isolate the systematic error induced by 
the pickup, for these simulations the only input signal is 
the Earth's magnetic field; the CMB and dust sky sig- 
nals are not included. These magnetic field time-ordered 
data are the inputs to the map-maker. In this case there 
is no differencing of the final map with an input map; the 
output stray field map is treated as the residual map. 

Figure |8] shows the resulting map for Q polarization 
signal induced by the Earth's magnetic field. The U map 
has similar amplitude. Figure [9] illustrates this magnetic 
response in multipole space. The power spectrum resid- 
uals, Ri, are less than 3% of the CMB i?-mode spectrum 
for r = 0.03 across all multipoles, and so we conclude 
that with the improved focal-plane shielding, the mag- 
netic pickup will not affect Spider's science goals. 

5.4. Mis-matched beams 

In this section we examine the impact of non-idealities 
in the beams produced by the Spider optical system. 
Ideally, the instrumental beams would be azimuthally 
symmetric and would be independent of the antenna po- 
larization. We consider several departures from this ideal 
situation, motivated by reported characterizations of the 
Bicep2 instrument (which shares a common optical and 
antenna design with Spider) and early measurements 
of the Spide r beams. In parti cular we use the results 
presented in lAikin et"al] ()2010f ) to guide us. Of most 
concern are mismatches between the beams seen by the 
two bolometers in a spatial pixel, as these have the po- 
tential to strongly couple total intensity on the sky into 
the observed polarization. In this work we consider three 
such effects that have proven to describe accurately the 
main departures from ideality in the measured beams: a 
differential pointing error of 1.5arcmin, a differential el- 
lipticity of 0.012 and a differential beam width of 0.5%. 
The method we employ to include these effects in our 
simulations is designed to be sufficiently flexible to be 
applied to a more detailed description of the beams, such 
as will be available when full beam measurements of the 
Spider instruments are carried out in the near future. 

Our baseline simulation method assumes that the in- 
strumental beams are ideal, i.e. that each bolometer 
couples to the sky via the same azimuthally symmet- 
ric Gaussian kernel. This azimuthal symmetry ensures 
that the computationally-intensive convolution operation 
need only be performed once for each sky pixel. For a 
general kernel this is not the case: the convolution is 
a function of the instrument orientation, which is con- 
stantly changing throughout the scan. Performing a 
convolution for each sample in the time-ordered data is 
not possible in a reasonable time. Therefore, to investi- 
gate the impact of non-azimuthally symmetric beams, we 
model them as the sum of a small number of Gaussians, 
each offset relative to the nominal pointing direction and 
with differing widths and amplitudes, 

N 

B{x) = ^ a„G(a;; a;„,CT„). (13) 

ri=l 

Here B(x) is the beam as a function of position on the 
sky when the instrument is in some fiducial orientation, 
G{x;xQ,a) is a unit-normalized Gaussian centered on 
a;o with width cr, and the amplitudes, a„, sum to unity. 
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Fig. 8. — Residual Q map for a simulation including the signal induced by the Earth's magnetic field in equatorial coordinates. The 
black lines show the Galactic cut used when calculating power spectra. The magnetic field for each Spider pointing was determined using 
the World Magnetic Model code and data provided by the National Geophysical Data Center. The resulting time-ordered magnetic field 
data were used to construct the observed Q map. 
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Fig. 9. — Residual spectrum, Ri, for a simulation of the signal induced by the Earth's magnetic field (green, long-dashed). The spectrum 
for our nominal simulation (red) and the CMB B-mode spectrum for r = 0.03 (black) arc shown for comparison. 



As well as tracking the location of the nominal pointing 
direction on the sky, the pointing generator has been 
modified to also track the location of the center of each 
sub-beam, so that the full beam can be reconstructed at 
all times in the appropriate location and orientation. 

To generate the bolometer outputs for this model, 
the contribution from each (azimuthally symmetric) sub- 
beam, d„, is first calculated using equation (jB]) or equa- 
tion p2|) . along with the sky-model convolved with a 
Gaussian of width (t„. As the convolution operator is 
linear, the final bolometer output is then given by 



N 



(14) 



This method only requires a factor of N more convolution 



operations (or less if some of the cr„ are equal to each 
other), as well as a factor of N more operations in the 
pointing and time-ordered data generation. Note that, 
in keeping with the general approach of this paper, no 
attempt is made to correct for these beam errors during 
the map-making and power spectrum analysis. 

So far, we have described a method to include non- 
azimuthally symmetric beams in our simulations. We 
also wish to consider the impact of differences in the 
beams seen by the two bolometers in each spatial pixel. 
This can be achieved by using different parameter values 
{N, (Tn and a;„) for the two beams, Bi{x) and B2{x). We 
can also allow the beams to vary across the focal plane 
by using different beam models for different focal-plane 
locations. 

We must be careful with the applicability of this model. 
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It is important that we only consider two different an- 
tenna polarization orientations across our focal plane. 
For other polarization orientations (e.g. at ±45 degrees), 
the instrument beams are functions of Bi{x) and B2{x) 
in a manner which does not generally conform to equa- 
tion (|13|) . Similarly, this model does not include (and 
cannot in general be extended to include) cross-couplings 
between the two polarization states induced by the opti- 
cal system. 

We perform simulations for each of the three beam 
mismatch errors we consider. In each case, we model 
the beams to hav e errors typical of those reported in 
lAikin et al.l ()2010D . This represents the worst case for the 
likely performance of the Spider instruments, as further 
research and development is expected to improve the fi- 
delity of the beams produced. Note that, in this section 
only, we slightly reduce the ideal beam fwhm used in the 
simulations to match those measured for Bicep2, from 
40 to 31 arcmin. 

We begin by considering differential beam pointing, i.e. 
differences in the beam centers seen by the two bolome- 
ters in a spatial pixel. Such errors couple gradients in 
the total intensity of the radiation incident on the opti- 
cal system into linear polarization in the focal plane. If 
such contaminat ion is not removed, the related tol erance 
limits are strict (|Hu et al.ll20"ollO'Dea et al.ll2007D . Spi- 
der was designed with such concerns in mind, however: 
the half-wave plate is on the sky side of all the optical 
components. The effect of this coupling on the bolometer 
outputs will thus be independent of the half-wave plate 
orientation, and so will be suppressed in the map domain 
to a degree dependent on the half-wave plate rotation 
st rategy. Note that unlike the pointing errors considered 
in iMacTavish et aD ()2008[ ) where the pointing properties 
were defined to be time-invariant on the sky, here they 
are defined to be time-invariant in the focal plane. 

Modeling the differential pointing is straightforward: 
we simply use one 'sub-beam' for each bolometer in a 
spatial pixel, with an appropriate offset between them. 
For this offset we use a typical value of 1.5 arcmin. The 
residual spectrum for a simulation using this model is 
shown in Figure 111! and the residual Q map is shown in 
Figure [T^] (the residual U map has similar amplitude) . 
On large scales, where the gravitational wave signal re- 
sides, the residuals are small, rising to around 20% of 
the CMB B-mode spectrum for r = 0.03 at I 100. 
On smaller scales the residuals are significant, becom- 
ing comparable to the CMB spectrum around / ~ 300. 
This does not compromise Spider's goals, however, as 
the gravitational lensing signal and instrumental noise 
are expected to dominate on these scales. Note that the 
residuals do not depend strongly on the orientation of 
the pointing error in the focal plane. Here we present 
the case where the pointing errors are parallel to the col- 
umn of detectors we simulate. Repeating the simulation 
with the pointing errors perpendicular to the column, the 
residuals are generally lower for I < 150, and similar on 
smaller scales. 

Comparing these low residuals to the strict tolerance 
limits found for differential pointing errors in the ab- 
sence of polarization modulation with a half-wave plate 
(|Hu et al.ll20"03l:IO'Dea et al.l[200l . we can conclude that 
our half- wave plate rotation strategy leads to strong sup- 
pression of these errors at the map-making stage. 



Next we consider differential beam ellipticity, which 
also couples the total intensity of the radiation inci- 
dent on the optical system into linear polarization in 
the focal plane, through a local quadrupolar pattern. 
Defining the ellipticity of a beam as e = (FWHM,„aj — 
FWHMniin)/(FWHM,„aj -I" FWHMmin) , where the subscripts 
denote the major and minor axes of the ellipse, we take 
a typical value of the difference in ellipticity between 
the beams for the two bolometers in a spatial pixel of 
61-62= 0.012. 

We model an elliptical beam using two sub-beams with 
identical widths, a = 31 arcmin, and centers displaced 
from the nominal pointing direction by 6/2 in opposite 
directions along the desired orientation of the major axis. 
We find that for S = 2 arcmin, e = 0.006. Figure [TOl 
shows the difference between this model and a true ellip- 
tical Gaussian with matching ellipticity; note that they 
agree to 0.1% everywhere. To produce the desired differ- 
ential ellipticities in our simulations, we use this model 
with the major axis parallel to the polarization direction 
for each antenna. The resulting residuals, shown in Fig- 
ure [TTl are negligible across the multipole range relevant 
to Spider. 

The relative orientation of the major axes of the el- 
liptical beams and the polarization of the antennas is 
also an important factor in determining the impact of a 
differential ellipticity error. In the absence of half-wave 
plate modulation, differential ellipticity contributes a sig- 
nal indistinguishable from a true polarization signal, i.e. 
one which cannot be mitigated through instrument rota- 
tion. The nature of this signal depends on the relative 
orientation of the major axes and the antenna polariza- 
tion: if the major axes of the beams are at 45 degrees 
to the corresponding antenna-polarization orientations 
this sig nal produces a -B-m ode pattern in the recovered 
maps (Sh imon et al.l 120081 ). Therefore, we also perform 
a simulation with the major axes of the elliptical beams 
at ±45 degrees to the antenna polarization orientations. 
This simulation produces residuals comparable to those 
shown in Figure 111! confirming that our half- wave plate 
rotation strategy is sufficient to mitigate the expected 
differential ellipticity in the optical system. 

Finally, we consider differential beam widths, i.e. dif- 
ferences in the widths of the beams for the b olometers in 
a spat ial pixel, taking a typical value from lAikin et al.l 
(20W) of 0.5%. To do so, we simply use one 'sub-beam' 
with the appropriate beam-width for each bolometer. We 
find that the residuals are negligible across all multipoles, 
with Ri remaining within 20% of its value in our baseline 
simulation; see Figure [TT] 

We have shown that the expected beam mismatch er- 
rors in the Spider optical system will not compromise 
the science goals of the experiment, even if no attempt 
is made to correct for them during data analysis, as 
the half-wave plate modulation works well in mitigating 
them at the map-making stage. Although we have not 
considered variations in the magnitude and orientation 
of these mismatches across the focal plane, such varia- 
tion is only likely to reduce the overall residuals, as the 
contributions from different spatial pixels combine with 
less coherence. 

6. CONCLUSIONS 
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Test of beam model 
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Fig. 10. — Map showing the difference between the elliptical-beam model we implement in our simulations, with ellipticity e = 0.006, 
and a true elliptical Gaussian. The beams are normalized such that the true elliptical Gaussian peaks at one. Note that the difference is 
small, peaking at around 0.1%. 
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Fig. 11. — Residual spectra, i?,;, for simulations with a differential pointing error of 1.5arcmin (magenta, dashed), a differential ellipticity 
error of 0.012 (green, long-dashed) and a differential beam width of 0.5% (blue, dotted). Note that the orientation of the pointing errors 
and ellipticities in the focal plane do not significantly affect the residuals. The spectrum for our nominal simulation (red), and the CMB 
B-mode spectrum for r = 0.03 (black) are shown for comparison. 

We have considered several potential sources of sys- 
tematic error in Spider and used simulations to assess 
their hkely impact on the experiment's main science goal, 
assuming that no attempt is made to correct the effects 
during data analysis. Unhke previous examinations of 
systematic errors, the goal of this work is not to set tol- 
erance limits to guide the design of the instrument, but 
instead to use measurements of particular instrumental 
non-idealities to assess their importance. Many of the 
systematics considered here have so far received little 
attention in the literature. We have also described a 
model of the polarized Galactic foreground emission Spi- 



der will observe, based on a three-dimensional model 
of the Galactic magnetic field and dust, and we have 
demonstrated that its presence will not compromise the 
performance of our baseline observation strategy. 

We have considered three main sources of systematic 
errors that have been characterized through measure- 
ments of Spider (or SPiDER-like) hardware; the non- 
ideal spectral response of the half-wave plate, spuri- 
ous signals generated by the motion of the instrument 
through the Earth's magnetic field, and mismatches in 
the instrumental beams generated by orthogonally po- 
larized antennas. 
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Fig. 12. — Residual Q map for a simulation with a differential pointing error of l.Sarcmin in equatorial coordinates. The black lines 
show the Galactic cut used when calculating power spectra.. 



We concluded that non-idealities in the spectral re- 
sponse of the half-wave plates result in errors in the re- 
covered B-mode power spectrum that may mildly bias 
our measurement of r if the gravitational wave ampli- 
tude is close to our target sensitivity of r = 0.03 or below. 
This systematic effect, left uncorrected, will not signif- 
icantly degrade the scientific results of the experiment, 
nor lead to a spurious detection of gravitational waves. 
However, it should be straightforward to remove this bias 
by introducing three parameters describing the half- wave 
plate non-idealities (which can be accurately measured 
pre-flight) into the pointing matrix during map-making. 

Our simulations of the stray magnetic field pickup in 
the TES detectors and SQUID amplifiers have demon- 
strated that the recent changes to t he focal plane design 
to improve the magnetic shielding (|Runvan et al.|[2010l) 
are sufficient to suppress the pickup to a negligible level. 
Using laboratory measurements of the focal plane mag- 
netic field response, including differential pickup within 
detector pairs, the simulated S-mode power spectrum 
residual is over an order of magnitude smaller than the 
CMB for r = 0.03 across the multipoles of interest. 

Finally, we have shown that the differential beam non- 
idealities measured in Spider and Bicep2, which has a 
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similar optical system, are not significant; the modula- 
tion of these errors introduced by the half-wave plate on 
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